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Abstract. We present a common framework to study decay and exchanges 
rates in a wide class of dynamical systems. Several applications, ranging form 
the metric theory of continuons fractions and the Shannon capacity of con- 
trained systems to the decay rate of metastable states, are given. 



1. Introduction 

In applications of the theory of dynamical systems to concrete situations it is 
often necessary to study rare events. Examples arc open systems with a small 
chance to escape and metastable states. Although much numerical work exists 
(e.g. see [(>, !)] and references therein) not many rigorous results are available. 
In principle one can try to apply perturbation theory but the existing theorems 
[13, 12] do not produce very sharp results. A similar situation occurs in the study 
of linear response theory. While perturbation theory applies to a wide class of 
smooth systems [21], this is no longer true when discontinuities are present in the 
system. In that case not only perturbation theory does not imply linear response, 
but in fact there are cases when linear response itself is violated, see [1, 3] and 
references therein. Since an open system is typically modeled by a hole in the 
system (that is by a region in which the dynamics stops once the trajectory enters 
it), the presence of discontinuities is inevitable. 

Accordingly, one could expect that the quasi-invariant measure (which describes 
the long time distribution of the trajectories conditioned to the event that they 
have not left the system, i.e. they have not entered the hole) and the escape rate 
(that measures the rate at which trajectories leave the system) depend in a very 
erratic (non-differentiable) way on the size and position of the hole. Yet, almost 
nothing is known about such situations. 

In the present paper we prove a general theorem providing a first order expansion 
of escape rates and exchange rates in terms of the strength of the rare effect that is 
investigated (whereby refining the results in [V-^] even-though limited to the present 
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setting). We derive from this theorem exphcit formulae for escape rates (both in 
one dimensional and two-dimensional cases) and for the exchange rate between two 
quasi-invariant sets (mctastability) . To our knowledge such formulae were known 
only for rather special cases and are completely knew in the generality presented 
here. 

Just to give an impression of the wide applicability of our main result, we list a 
few examples that are detailed in section 3. 

Let z be a periodic point of period p for the doubling map x ^—^ 2x mod 1 and 
consider intervals 1^ 5 z of length e. Denote the decay rate for the hole by A^; 
i.e. the Lebesgue measure of the set of points that are not trapped by during 
the first n iterations of the map decreases asymptotically like A". Then A^ has the 
following first order expansion at e = 

(1.1) A, = l-e-(l-2"f)+o(£) . 

Similar formulae can be obtained when the hole is a union of several intervals, and 
analogous results hold of course for a coin tossing process that is stopped once a 
pattern of heads and tails from a given finite collection is observed. 

We turn to continued fraction expansions. Let mfc,n be the Lebesgue measure 
of those points x G [0,1] whose continued fraction expansion up to the n-th digit 
does not contain a block of k consecutive I's. For each k, these numbers decrease 
asymptotically like some AJJ and, setting z ~ , we show 

(1.2) l„i_|i = £!0_ti!)!. 0.6504. 

k->oo z^*^ in 2 

The above are one dimensional examples. Along the same lines one can treat 
piecewise expanding maps in higher dimensions provided the invariant density is 
not too irregular in a neighborhood of the holes. This is not guaranteed working 
only with the usual multivariate BV-spaces, but rather by variants as developed by 
Blank [4] or Saussol [23]. In principle the present theory also applies to Anosov dif- 
feomorphisms if one can devise the proper functional space setting. Unfortunately, 
despite recent progress [S, 2], the available settings are still not adequate for the 
applications of the present results. Nevertheless, it is conceivable that in the near 
future this result could be applied e.g. to billiards. 

A related but different question occurs if a system has two ergodic mixing com- 
ponents that share a common part of their boundary in phase space. In that case 
small perturbations can cause rare "jumps" over the boundary giving thus rise to 
quasistationary (also called metastable or nearly invariant) behavior. As a result 
the double eigenvalue 1 corresponding to the two original mixing components splits 
into a single eigenvalue 1 and another real eigenvalue Ag close to 1 which char- 
acterizes the rate of exchange between the two components, see e.g. [7, 20]. In 
subsection 3.2 we consider piecewise expanding ID-maps T on the interval [0, 1] 
with two mixing components Ii , I2 having a fixed point z in common. The Markov 
process obtained by adding to the dynamics, at each time n, independent identically 
distributed random noise eZ„ shows quasistationary behavior with 

(1.3) Mm'—^^^fi- n\w\] + ^nz,] 



6^0 e 2 V T'{z) J " " 2 

where a = {2m{h))-\ /3 = (2m(/2))"i and W := Er=o[^'(^)]""^"+i- 
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The paper is organized as foUows: in the next section we describe the general 
setting and state our main Theorem 2.1 whose proof is postponed to section 6. In 
section 3 we apply the theorem to decay and exchange rates of piecewise expand- 
ing interval maps and illustrate its applicability with some specific examples: the 
doubling map, the Gauss map and the generalized cusp map. The decay rate of a 
two-dimensional example is studied in section 4. More precisely, we study the rate 
at which trajectories of two coupled ID-maps synchronise up to some difference e in 
the limit £ ^ 0. Finally, in section 5, we indicate relations among between formula 
and approaches to metastability in molecular dynamics [1!)], in oceanic structures 
[11], and with the Shannon capacity of constrained binary codes [1 ')]. 

2. An ABSTRACT PERTURBATION RESULT 

Let (y, ||.||) be a real or complex normcd vector space with dual ||.||) . Con- 
sider a family : V V [e ^ E) oi uniformly bounded linear operators where 
i? C K is a closed set of parameters with e = as an accumulation point. We make 
the following assumptions on the operators P^: there are g C, <y9e G i^e G 
and linear operators Q^, -.V -^V such that 

(Al) A^lPe = 93. Z^e + Qe, 

(A2) P,(^,) = Ae(^e, VePe = QeiVe) = 0, iy,Q, = 0, 

OO 

(A3) VsupllQ^'il =: Ci < oo, 

(The summability condition in (A3) can only be satisfied if the operators Pe have a 
uniform spectral gap. See Remark 2.2 for a weakening of this requirement.) Observe 
that assumptions (Al) and (A2) imply Ve{fe) = 1 for all e. As our ultimate goal 
is to prove a perturbation result for small e, it is natural to relate the "size" of (^^ 
to that of LpQ by a further assumption: 

(A4) vai^e) = 1 and sup \\ipe\\ —■ C2 < 00. 

Finally we denote 

(2.1) A, := uoiiPo ~ PeK^o)) 

and we make the following assumptions to control the size of the perturbation: 
there is C3 > such that 

(A5) r;e:= 111^0(^0-^5)11^0 as e^O, 

(A6) r;,.||(Po-Pe)(^o)il <C3|A,|. 

Here 77^ denotes the norm of the linear functional i'o{Po — P^) V ^ M.. 
The basic identity is 

(2.2) Ao - Ae = XoM^e) - M^ei^e)) = M{Po - Pe){^e))- 

In view of assumptions (A4) and (A5) this implies 

(2.3) |Ao - A,| < C2Ve, 

in particular, lim^^o A^ = Ao- The main result of this section is the following more 
accurate approximation for Ao — A^. 

Theorem 2.1. Assume (Al)~(A6). 
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a) There is eo > such that = Aq if s < eq and Ag = 0. 

b) If 7^ for all sufficiently small s Cz E and if 

(A7) q, hm q,, := hm 

exists for each integer k > 0, then 

^ k=0 

Remark 2.2. In section 6 we prove this theorem under slightly weaker hypothesis 
that may be applicable also in non-uniformly hyperbolic situations. Namely, we relax 
the summability condition X]J^°=o ^'^PeS-E HQ" II < °° from (A3) in the following way: 
there is a second norm ||.||* > |M| onV such that 

oo 

(A3*) Vsuplig^r ==:Ci<oo 

where ||Q"||* := suplllQ"?/;!! : HV'H* ^ !}• We have to compensate this by slightly 
stronger assumptions on ipo, namely \\^q\\* < C2 < 00 and 

(A6*) Ve-\\iPo~Pe)MU<C3\A,\. 



3. Applications to piecewise expanding interval maps 

Assume that T : [0, 1] [0, 1] is piecewise monotone with (possibly countably 
many) continuously diffcrcntiable branches. (This means that each branch is con- 
tinuously differentiable in the interior of its domain so that the derivative even of 
a single branch may be imboimdcd.) Define g : M — > M by g(x) = l/|T'(a;)| if x 
is in the interior of one of the monotonicity intervals of T and g{x) — otherwise, 
and assume that \\g\\oo < 1 smd that g is of bounded variation. Let BV be the 
space of real- valued functions of bounded variation on [0, 1]. Rychlik [22] showed 
that the Pcrron-Frobcnius operator P of T acting on Lebesgue equivalence classes 
of functions from BV is quasi-compact. (As BV- functions have at most countably 
many discontinuities, two BV-functions in the same Lebesgue equivalence class 
have the same discontinuities and differ at most by their values at these countably 
many points, and we will not distinguish henceforth between BV-functions and 
their Lebesgue equivalence classes.) 

If T is mixing this implies that Pq = P satisfies (A1)-(A3) for e = with 

= Lebesgue measure), Aq = 1 and < & BV. 

The essential observation behind this is that for e = a Lasota-Yorke type 
inequality is satisfied: there are constants r G (0,1) and R> such that for 
e = 0, aU n e N and aU / e BV, 

(3.1) ||P;7ll<i?(r"||/||+ I I/I dm) 

where ||/|| is the variation of the extension of / to the whole real line by setting 
f{x) = if X ^ [0, 1]. For the applications in this section we will assume that this 
inequality holds not only for e = but, with uniform constants r and R, for all 
£ € E. This is mostly the case when P^ is a small dynamical perturbation of Pq - 
however there are exceptions, see [13] for a more precise discussion and references. 
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3.1. Decay rates. We suppose that T is mixing. Let (V, ||.||) be the space BV, 
let E = [0, El], and consider a family {Ie)eiEE of compact subintervals of [0, 1] such 
that C 4, if e < e'. Define the operators by P^if) = -P(/l[o,i]\/J- ^ 
m{Ig-^) is sufficiently small, the perturbation results from [13] apply provided (3.1) 
holds, see [18, section 7]. In particular, (Al)-(A4) are satisfied for e £ E. We have 
= i>o{P{lj^\jgipo)) ~ ^q(I^\Iq) where fiQ is the probability measure with density 
ipo w.r.t. vq. (fiQ is indeed the equilibrium state for logg on the "non-trapped" set 

o 

Xnt ■= {x e [0, 1] : T"x ^ /qVti > 0}. (pe is also the conditionally invariant density 
for the "hole" I^, see e.g. [L"^].) 

We need to check assumptions (A5) and (A6). First note that 



(3.2) 77^= sup |t/o(^o(V'l/,\/o))l = l-^ol sup 

||i/j||<i ll'i/'ll<i 



< \Xo\Mh\Io)^ 



In particular, \i^o{Po~Pe){fo)\ = fodi^o. As \\{Po-P^){(po)\\ < O{\\ipoii,\io 

assumptions (A5) and (A6) will be satisfied if vo{Ie \Io) ^0 when e — > and if 

(3.3) ||y^olj,\Joll < const \ / (fiodi^o- 

^OVe \ ^0) Jl,\Io 

This condition (as well as conditions (Al)-(A4) discussed above) can be checked 
easily in most cases of interest. It is always satisfied if inf (/JqI/ej > 0. 

3.1.1. Holes shrinking to a point. We specialism to the case where Iq ~ {2} for 
some z £ [0, 1] so that Pq is indeed the Perron-Frobenius operator P for T and 
Aq = 1, and we assume for simplicity that T and also the invariant density ipo are 
continuous at z. We consider with length e, so TO(/e \ Iq) = and we assume 
that > 0. Here are a few examples: 

The doubling map: T{x) = 2x mod 1 with </Jo(a;) = 1. 
The Gauss map: T{x) = ^ modi with ifo{x) = j^j^ 
The generalized cusp map: T^{x) = 1 — \2x — 1|'^ for some 7 G (|, 1]. As 
\T'^{x)\ ~ ^2x-i\'^-i — ^^^^ map is a uniformly expanding map with 
two full branches. The weight function g{x) = |T^(a;)|^^ has two monotone 
bounded branches and is clearly of bounded variation. (Observe that Ti 
is just the tent map. T1/2 is known as the cusp map. It has a; = as a 
neutral fixed point and is not covered by the present setting.) The invariant 
density foix) of behaves like ^£2ilZHl(i _ x)^~^ near a; = 1, so it has a 
zero at a; = 1 if 7 < 1.^ 

In all three examples, A^ = f-oih) > 0. In the first two examples, condition (3.3) 
is clearly satisfied because inf (po\i^ > if e is sufficiently small. For the generalized 
cusp map the same is true if z ^ 1. In case z = 1, ||(/3ol/e|| — 2const^e~~^ and 
J J ipo dm = const^ Ji^^i^ — x)~^^dx = const-^ 7 e~ so that (3.3) is satisfied as well. 



^Here is a sketch of the argument: As the map has full branches, the invariant density (po = 
limn_,oo ^0*1 is continuous. Also ipo(^) > 0, because otherwise 1^0(2:) = for all x g (J^ '^-^"{f }' 
and this set is dense in [0, 1]. Therefore, for x close to 1, 

ipo{x) = Po^poix) ~ (1 ~ x)y 
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Let 

Uk,e ■■= T-\[{), 1] \ 4) n • • • n r-'=([o, i] \ h) n r-^^'+^h,. 

As qk^e = Mo(^£ n Uk,e)/fJ.o{Is), we find: 

If z is not periodic: then Uk,e = for sufficiently small e so that qk = 
liniE^o Qk,e = for all k. Therefore, 

(3.4) lim — - — ^ ~ 1, in particular lim — ; — f ~ foiz) ■ 

If z is periodic with period p: then Uk.e = for sufficiently small e except 
if A; = p — 1 so that 

(3.5) liml^ = l-lim^"^^-"jr^-)^l ' 



.0^0 (4) Mo (4) 

(3.6) in particular lim — r— ^ = <^o(^) f 1 



om(/,) ^-^V l(TP)'(z)|, 
Formulas (3.4) and (3.6) imply that the function e Ae is differentiable at e = 0. 
Note, however, that in general it is non-diffcrentiable at other values of e, see section 
3.1.2 below. 

We look more explicitly at the above three examples. Recall that Iq ~ {z}. 

The doubling map: lim^^o ^ 1 - 2-p if TP{z) = z. This is (1.1). 

The Gauss map: a) Consider z = and e ^ E := {i, i, i, . . . } U {0}. Let 
Je := [0, e]. In terms of the continued fraction algorithm this means that 
the expansion stops as soon as a digit > is generated. In this case it is 
easy to see that 

n Uk,e) < n T-C^+i)/,) = P^iPoliJdm 

= 0(e/io(4)) 



so that qk^e — 0{e) and hence qk = for all k. Hence lim^ ~ ^ — i 



ln2 ■ 



b) For the same map we consider z = which is the rightmost fixed 
point of T . We have T'(z) = — z^^. As —z and z~^ are the two zeros of 
cc^ — a; — 1, it is obvious that 1 — m/, = z. Hence, for intervals 1= of length 

£ around z we have limg^^n ~ 7^ = t-^ ■ 

c) Denote /(x) = ^ — 1. Then / is the rightmost branch of T, and the 
interval around z which is mapped by onto (0, 1) has endpoints 

and /~'^(0). Denote the length of this interval by and the interval itself 
by /e^. As —z and z"^ are the eigenvalues of the the coefficient matrix of 
a calculation shows that = z'^^^'^{\^ z'^)'^{\ + 0{z^^^'^))? In terms 
of the continued fraction algorithm the hole 1^^. means that an expansion 



^Denote the coefficient matrix ("J) of /-^ by M and let =: ("'=^''). Then f-''{xi) - 

f-H^o) = <ietiM'^) ,.,.,+Tj,::.,+,,, so that = - f->=7o)\ = J^^^^r 

Cfc+i = dfc and d^+i = (c^ we have = , . ^ . , and d^+i = +dfe_i. With dg = di = 1 

I 

this yields dk = (^-(^+2) + {-z)*"). Hence 
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stops at time n + k as soon as at least k consecutive digits 1 are generated. 
So it is natural to rewrite the limit from b) as in formula (1.2), namely 

3.7) lim ^ = \ ~ 0.6504 

The generalized cusp map: We focus on z ~ 1, where the invariant den- 
sity (pf) vanishes, and consider holes = [1 — s,l]. Then, for e close 
to 0, we have ^o(^e) ~ fo{i/2) J^^^(^i _ x)~^^ dx ~ (po{^) e~ so that 

3.1.2. Holes shrinking to a nontrivial hole. We assume now that /q is a nontrivial 
interval of some fixed length ^ > and the intervals /g 3 /q have length I + e. To 
simplify the discussion we assume more specifically that = [a — £, a + ^] where 
a S [0, 1] is a continuity point of T and also of (/Sq- Assume furthermore that a is 
not periodic for T (the periodic case can be dealt with analogously). Now, as /zq is 
supported by the non-trapped set Xnt, we have in particular ^j.o{Io) = and hence 
= fio{I^). It follows from Theorem 2.1 that either fJ.o{Ie) = and hence = Aq 
for all sufficiently small e, or lime^o ^°7/^ ^ ^- observe that fiQ is of fractal 



nature, so typically /io(^e) depends on e in a devil's staircase manner. Hence cither 
Ae = Ao for small e (which happens if a itself is trapped) , or linie^o ^^'(^j ) — and 
e ^ Ae is not differentiable at £ = 0. 

3.2. Exchange rates. We suppose that T has two ergodic components and that 
its restriction to each of these components is mixing. So the eigenvalue 1 of P has 
(geometric) multiplicity 2 and the rest of its spectrum is contained in a disk of radius 
smaller than some 7 G (0, 1). Let (ng)^^^ be a family of Markov operators close 
to the identity with Hq = I, and denote := P o n^. (One could as well consider 
Ilg o P since that operator has the same eigenvalues as P^.) Under rather weak 
regularity assumptions on the He, the spectral perturbation results from [13] apply 
again. This is true, for example, if the He are convolutions with smooth densities 
ke{x) = e~^k{e~^x) (modeling random perturbations) or if they are conditional 
expectations w.r.t. m and a finite partition into intervals of length s (modeling 
Ulam's discretization scheme.) As the Pg are also Markov operators, this means 
that 1 is an isolated eigenvalue of each P^. If it has multiplicity 2 there is nothing 
more to say about it. If it is a simple eigenvalue, however, then there is a second 
simple eigenvalue A^ close to 1 to which we will apply Theorem 2.1. 
Let {V, \\.\\) be the space 

BVo := {/ e BV : m(/) = 0}. 

BVo is invariant under all P^ , and the previous discussion implies that assumptions 
(Al)-(A4) are satisfied. More precisely, Aq = 1, and there is an increasing function 
X ■ [0,1] { — 1,1} such that x o T = x, \tpo\ — Xfo is an invariant density 
for P = Pq, and i^o = X'^^i so that ^0 '■= Vo^'o = I'Pol"^ is an invariant probability 
measure that gives equal mass to both ergodic components of T. Let /i = {x = — 1} 
and /2 = {x = 1} be the two invariant components of T . 
In order to apply Theorem 2.1 let 

Pe{x) ■.= \{l-{x{x)-{Iilx){x))) 
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where 11* is the dual of He with respect to Lebesgue measure on [0, 1]: It is easy to 
see that pe (x) is the probability, that the Markovian dynamics He move the system 
from the state x (that belongs to one of the two invariant components of T) to 
some state in the other component. 

Wc have to check assumptions (A5) and (A6). Routine calculations show that 
t'o(-Po ~ Pe){'4') — 2to(x ■ ps ■ i.') for each ijj G V. This implies 

(3.8) A, ^ vo{Po - Pe){vo) - 2m(|^o| ■ Pe) = 2/^o(Ps) 

and vq{Pq — Pe){i^) = '^'^ixPs'^) < |Pe|(i?Ti so that < 2m{ps). Therefore 

we require that the average probability m{p^) to change the invariant component 
under the action of II^ tends to as e — > and that 

(3.9) 11(1 - He) (y^o) II < const / \^o\Pe dm . 

m{pej J 

In the following wc will assume 

M\ipo\ > 0. 

This trivially implies (3.9) although the latter can be verified in many other cases. 
It remains to check assumption (A7). Observing (3.8) and 

'^0 {{Po - Pe)Pe'{P0 - PeKvo)) - ™((x - H^x) ' ^'(^0 " Pe)M) 

= 2m{p,x-PePoi^o-Tl,^o)) 
we get the following expression for the qk,£- 

(3.11) qk,e - ,1 \ , m (p, X ■ PePoi^o - n,v5o)) 

m(|(^o|Pe) 

The evaluation of the limit as e ^ depends strongly on the details of the map T 
and of the perturbation. We therefore make some further simplifying assumptions: 

• The Ilg are local perturbations, i.e., for each x, HgSx is supported in a 
Ce-neighborhood of x. 

• T is continuous. 

As the restrictions of T to its two ergodic components are mixing, the continuity 
of T implies that the non-wandering part of these components are just two single 
intervals li C h. If these intervals do not have a common end point, then = 
2Mo(Pe) = for small e so that = 1 for such e by our main theorem. Otherwise 
Ii and I2 have a common endpoint z. Since two interval can have at most one 
common endpoint and since the map is continuous, it follows by the invariance 
of /i,/2 that z is a fixed point. In this case, Pe{x) = unless x belongs to the 
Ce- neighbor hood of z. As an example let us consider the special (but still rather 
general) class of examples characterized by the following properties 

. ITUT; = [0,1].^ 

• Assume Ili,f{x) — K^{y,x)f{y)dy where is a positive kernel such 

that, for all y G [0,1], Ki.{y,x)dx = 1. In order to satisfy assumptions 
(Al)-(A4) one should suppose that the kernels are bistochastic or close to 
convolution kernels in the sense of [5, Corollary 3.20]. 



If the wandering part is present, the final result still holds with Ii substituted by UngNT"/; 
in (3.13). 
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• There exists a > such that i^e(y, x) = e^^ K{e^^{x—y)) provided \z~y\ < 
a. Here K is a smooth probabihty density supported in [—1, 1].^ 

• ipo is continuous in each crgodic component and T is diffcrentiablc at z. 

Note that, since \ipo\ must give the same weight to the two ergodic components 
it wiU, in general, be discontinuous at z. Let a, /3 be the left and right limit 
respectively. Then, introducing coordinates x ~ z + e( and setting 9{y) ~ signy, 
we have (po{z + eC) = ^^9{() + + o{l) , uniformly for in a compact set. Next, 
for e small enough, 

Peiz + e0^l f K{y - C)[l - 0{y)e{(:)]dy p(C). 

Note that p > and p(C) = if |C| > 1- Accordingly, 

Mo(Pe) = ■m{\ipQ\pe) = e^-^m{ep) + £ ^ ^ " m(p) + 0(e) =: eV + o{e). 

In addition, for each function ^'^ such that '^^{z + eC,) = ip{C) + for some fixed 
compact support function ip, holds, in the hmit e 0, (ne*e)(z + eC) — /r^(C'~ 
C)V'(C'X + 0(1) and {Po'fe){z + eO = A- 7/;(AC) + o(l) where A := > 0. 

Hence, 

(P,*e)(^ + eC) = A / K{CA~CMC)dC + o{l) 
=: (/C^O(C)+o(l). 
The above setting applies to 'I'g ^ Po{I — Il^)Lpo, namely 

[Fo(I - n,)H(^ + eO = ^[A(?(C) - 0(C)] + 0(1). 
Indeed, IC0{C) = A0(C) for |C| > A^^, hence KO — K.9 is compactly supported. Thus 

* ~ 2f ■ 

Since the operator K, has L°° norm smaller than |A|, the latter equality implies 

lini ^ 2r - f;(20p, ^"(AI - m- 

^ n=0 

On the other hand, a direct computation (observing the fact that d{y) = 9{Ay)) 
shows that 2ep = (I - /C*)6', so 



lim^"'^-^2r-^ 
e^o e 2 



, (AI - IC)9) - lim ((/C*)"6i, (AI - IC)9) 



^One can consider the more general case K^{y,x) = e ^K(y,e ^{x — y)), for some smooth 
function K. The final formula then holds with K{-) replaced by K{z, ■). 
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Note that (/C*)"^ converges pointwise to a function 600 such that 9 — Boo is supported 
in the interval [— (1 — A)^^, (1 — A)^^], see (3.15). In particular, lC*6oo ~ 9oo- Then, 

lim = 2r - '--^{0 - (AI - IC)9) 

- +0 e 2 

(3.12) 

[{2ep, 9) + {A- l)m(l 

^ ^ mil ^ 99^) + (^—^m{29p) 



2T - 


P 


+ a 




2 


2T- 


P 


+ a 




2 


2T - 




+ a 




2 


(3 + 


a 1 




2 





T'{z), 

To make the formula more explicit and transparent let us make some further 
remarks. The original dynamical system has a natural invariant measure defined 
by /i = lim„^oo ^ Sl-=o ^o'-'^- assumptions h is continuous at z and {/i, xh} 

form a basis for the eigenspace of the eigenvalue one of the operator Pq- Thus (po = 
+ ^^^h. Remember that a, (3 are chosen so that J^^ |(^o| = //^ Iv'ol c'™ = 
^. Hence ^ = a Jj^ h dm ~ a lim„^oo ^ X]fc=o /o^ ^^T^ ° "^"^ ~ am{Ii) so that 

(3.13) a = (2m(/i))"^ and, analogously, [3 = (2m(/2))"^ . 

To describe the meaning of the two factors involving 9 and 9oo , let Z he a random 
variable whose distribution has probability density K. Then 

(3.14) mi29p) = ~2E[Z] . 

Next let Zi, Z2, . . . be independent copies of Z. The kernel /C* describes a Markov 
process 

X„ = A-i(X„_i + Z„) - • • • - A-" 1^X0 + A'^-i^,^ . 

The asymptotic behavior of the process (Ar„) is determined by the random variable 
W ■= 'ET=i^''~^^k- Indeed, let Xq = C- Then A:„ -> +00 if > -C and 
Xn —00 if ly < (As the Zk have density, W ^ has probability 0.) As 
((/C*)"6')(C) is the conditional expectation of 9{Xn) given Xq = C, it follows readily 
that 

9^(0 = P(X„ -> +CX)|Xo = C) - P(^n ^ -^1^0 = C) 

= 1 -2¥{W < ~C) ■ 



(3.15) 
Hence 

(3.16) to(1 - 99oo) = 2 [F{W > C) + ^-W > C)] dC = 2E[|M^|] . 

Jo 

Note that {1 - l/T {z))m{l -~ 99^) = 2(1 - A)E[|iy|] > 2(1 - A)|E[VF]| = 2|E[Z]| = 
\m{29p)\, so the r.h.s of (3.12) is clearly positive. 

We finish this section with a comment on the term ^exchange rate. Let Af := 
{ifie > 0}, A- := {ip, < 0}, and Pe := 1^+ ■ P:1^- + 1^- ■ Pn^f- Pe{x) is the 
probability to exchange the sets under the action of P* . Now Proposition 5.7 
from [7] can be rephrased in our setting as 1 — = 2 J Pe\'^e \ dm,^ so it is nearly 
twice the "stationary exchange rate" J Pehg dm where /i^ = Peh^ is the unique 
invariant probability density of the perturbed system. If all Ai^ are identical (e.g. 



^The operator P in [1] corresponds to our and the signed measure u to our ip^m. 
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under suitable symmetry assumptions on the system as in [i , Corollary 5.9]), then 
coincides with from above. 

4. An APPLICATION TO TWO COUPLED INTERVAL MAPS 

Let T : [0, 1] — > [0, 1] be a mixing piecewise expanding map as in section 3. To 
simplify the discussion we assume that 7 := inf |r'| > 4. Let M := [0,1]^ and 
define, for (5 G [0, ^ — i), the two-dimensional coupled map 

f:M^ Af, f{x, y) = ((1 - 5)T{x) + 5T{y), (1 - 5)T{y) + 5T{x)) . 

It is uniformly piecewise expanding with minimal expansion strictly larger than 2 
in the sense that 

As discussed in great detail in [14] there is 5i G (0, ^ — ^] such that, for 5 G [0, 5i], 
T is mixing in the sense that its Perron-Frobenius operator P : BV(A/) BV(Af) 
has a unique invariant probability density h and a spectral gap. Here BV(A'f) is 
the space of functions of bounded variation on that vanish outside A/. 

For e £ E := [0, ei] let := {{x,y) & M : \x — y\ < e}. If we interpret 
Se as a hole in the phase space M, this means that we stop a trajectory as soon 
as the two components have synchronized up to a difference of at most e. The 
corresponding Perron-Frobenius operator : BV(A/) BV(A/) is defined by 
Peiip) = Pilp • 1a/\Sj)- Denote the (two-dimensional) variation of a function ip G 
BV(Af) by Var('0). It is easy to check that Var(Pe?/') < 2Var(F'0) so that the 
family of operators P^ satisfies a uniform Lasota-Yorke inequality. (Observe that 
we made the generous assumption 7 > 4 and consult [14].) In view of the spectral 
stability results of [lo], assumptions (Al) - (A4) are satisfied with vq = m (the 
Lebesgue measure on Af), tpo = h, = hm and Aq = 1. 

We turn to assumptions (A5) and (A6). Observe first that 

(4.1) MPo - Pe)W = mi^pisj < Ce Var(V') . 

(The constant C depends on the details of the definition of the variation.) So 
in particular rj^ < Ce and (A5) is satisfied. As h is of bounded variation, we 
may assume that it is regularized along the diagonal of M in the sense that for ID- 
Lebesgue-almost every x the value h[x,x) is the average of the limits oih{x—u, x-\-u) 
and h{x + u,x — u) as it \ 0. In view of (4.1) we therefore conclude 

(4.2) lim(2e)"^Ae = lim(2e)"^ / hdm^f h{x,x)dx. 
-^-o -^-"0 Js^ Jo 

As Var(/il5^) < 2 Var(/i), we conclude that (A6) is satisfied if h{x, x) dx > 0. It 
remains to evaluate the qk- As in section 3 let 

f/fc^e ■.^f-\M\Se)n---nf-''{M\Se)nf-'^''+^^Se . 

Then ^ = /^o(5'e H Uk,e) / l^o{Se) and, since the diagonal of M is invariant under 
T, we find 

I - 1 
qa = lim 90, e = -TiT', TT / ^(^' ^) ~ ^<;\\rrii^\\ 



jlh{x,x)dxJo ' {l~mT'{x)\ 
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and Qk = linij^o 9fc,e = for all fc > 1- So finally, 



1-A. . 1 



(4-3) 1"%^ = X ^^"'"^ V (l-2^)|nx)| ,' 



5. Related results 

5.1. Metastable states in molecular dynamics and oceanic structures. 

Phase space methods to characterize biomolecular conformations as metastable 
states are used in molecular dynamics (see e.g. [19] and references cited there). Very 
roughly, if the Markov operator describes the discrete time evolution of such a 
system in a fixed time scale and ii X = D1UD2 is a decomposition (up to null sets) of 
the underlying phase space, then the metastability measure of this decomposition is 

defined as meta(Di, D2) = ^ ^^(-Di)" V^,^ -Pel_Did/^e + ^£(-^2)" ■ 
Theorem 1 of [19] relates meta(I?i, D2) to the second eigenvalue of in a way very 
similar to formula (3.12). Reference [10] is an up-to-date review of the phase space 
decomposition approach to metastability in general flow dynamical systems, and 
[11] is an application of these ideas to the detection of coherent oceanic structures. 
In the framework of weakly coupled rapidly mixing Markov chains, reference [20] 
also relates the second largest eigenvalue of a system to the exchange probabilities 
between its components. 



5.2. Shannon capacity of constrained systems of binary sequences. In in- 
formation theory, the topological entropy of subshifts of {0, 1}^ that arc determined 
by a (short) list L,„ = {b[™\ . . . , Bp™'') of distinct blocks of length m which are 
not allowed to occur [17] is called the Shannon capacity of the system. It is closely 
related to the rate of periodic prefix-synchronized (PPS) codes with markers B^"^\ 
i = 1, . . . ,p (see e.g. [15]). For each sequence Li,L2, ... of such lists with fixed 
length p there are a subsequence {nij) and zi, . . . , Zp G {0, 1}^ such that the B^^^'^ 
converge to Zi as j — > 00. We will assume without loss that the full sequences 
(B'f^^)m converge. As the full two-shift is isomorphic (for each invariant measure 
of positive entropy) to the doubling map. the shift constrained by the forbidden 
blocks in is isomorphic to the doubling map T with "hole" Im being the union of 
those monotonicity intervals of T™ labeled by the words in L„i . Hence the topolog- 
ical entropy h{h,„i) of this shift equals log(2A„i) where A,„ is the leading eigenvalue 
of the Perron- Frobenius operator of T with hole /„, compare subsection 3.1. 

If the limit points Zi belong to b\"^^ for all i and m, one can analyses the 
situation just as in subsection 3.1. Some elementary reasoning yields the following: 
for i = 1, . . . ,p let = min{j > 1 : Zi G {zi, . . . , Zp}} with the convention that 
— +00 if no such j exists. Then 

hm l^gl_M!^^y(i_2-.W). 

m — >oo 2 — ^ 

i=l 

This is minimal when all = 1, e.g. if the Zi form just one periodic orbit. In 
that case h{hm) = log 2 — p2^^™'^^^ + o(2~™), which supports the conjecture on 
the precise values of h{hm) for p and m > fc 1 stated in [15]. 
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6. Proof of the main theorem 

As announced in Remark 2.2 we prove Theorem 2.1 under the weaker summa- 
bihty assumption (A3*). The reader who does not want to follow this slight gener- 
alization of the argument may just neglect all attached to the norms. We use 
the following notation: 



^ sup II g;^ I 



Lemma 6.1. There is a constant C > such that, for all e ^ E and all N > 0, 

a) |1 - Vei'Po)\ < C rje, 

h) IIQfy^oll <CKAr(||(Po-Pe)(<^o)||* + |Ao-Ae|). 

Proof, a) As (I - \-^Pe){\-^Pef{fo) = (I - K^Pe)Ql{fo) for all fc > 0, 
|l-z.,(^o)l= lim \ya{vo-{\-^P,r{vo))\ 

oo 

<^|z.o((i-Arip,)g,n^o))| 

oo 

= Y.h{i>^o'Po-K'Pe)Q';{^o))\ 

oo CX) 

< lAorv WQeirM* + iAo 'iiAe - Aoiikoii E WQeiryoW* 

fe=0 k=l 

= 0{Tje) + 0{Xo - A,) = 0(rk) 



where we used (A3*) and (2.3) for the last estimate, 
b) For each > wc have 



IIQf^oll <limsup||gf (^0 - (A7iPe)"(^o))|| +limsup|igf+>o|| 

n — >-oo n— i-oo 
oo 

<EllQe"(A7'-Ps)'(I-A7ip,)((^o)||+limsupKjv+„||<^olU 

OO 

<|Ao 'I E iWQ^^'iPo - Pe){^0)\\ + |Ao - Aelllgf +^+Vo| 



fc=0 



<|Ao 'I E WQ^^'W* (11(^0 - Ps)(^o)||* + |Ao - Aelll^olU) 



k=0 



=0(«;w)(|j(Po-Pe)('^o)||* + |Ao-Ae|) 



□ 
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Proof of Theorem 2.1. Observe first that by (2.2), for each n > 0, 

^ei^o) (Ao - K) 

=A, - ^ ((^^0 - P.)(A-ip,)^(I - X-'Pe){^o)) + 0{v,\\Q^^o\\) 

k=0 

Tl-1 

= A, - Ao 1 {{Po - Pe){Xe'Pe)''{Po - Ps)i^o)) 

71 

+ X^\Xo - Ae) {{Po - P,){Xj'P, f{^o)) 

fc=i 

+ C'(k„)(|A^| +r?e|Ao - Ael) (by Lemma 6.1b and (A6*)) 

=A, (l-Ao^xiArO 

n 

+ 0(%)|Ao - A,| J2 iWe{^0)\ye\\ + llQeVoll) + 0{k,,){\A,\ + 7J,\Xo - X,\) 
k=l 

where 

.... 1^0 {{Po - Ps)PHPo - Pe){V0)) 

^ ^ ^0iiP0-Pe){^0)) 

Observing Lemma 6.1a, (A4) and Lemma 6.1b, the error terms can be estimated 
by 0{ri^) n\Xo — A^l + C'(K„)|Ae| so that, in view of Lemma 6.1a, this yields, for 
each n > 0, 

(6.2) {l + 0{rj,)){Xo~X,){l + nO{v,)) = A, ( 1 - A^^^ |] A^'^gfc J +0{nn)\Ae\. 

If Ag = and is small, it follows that A^ = Aq. Otherwise we assumed in (A7) 
that (7fe = linij^o 'ik,e exists for each /c, and we conclude 

lhn^^^ = l-X:V''-^^''Z. + 0(«n) 
for each n > 0. From this the claim (2.4) follows in the limit n ^ oo. □ 
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